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Using the most recent cosmic ray spectra up to 2x10 eV, production 
spectra of high energy v's from cosmic ray interactions with interstellar 
gas and extragalactlc interactions of ultrahigh energy cosmic rays with 
3K universal background photons are presented and discussed. Estimates 
of the fluxes from cosmic diffuse sources and the "nearby" quasar 3C273 
are made using the generic relationship between secondary v's and Y's 
and using recent Y -ray satellite data. These Y -ray data provide impor- 
tant upper limits on cosmological v's. Quantitative estimates of the 
observability of high energy v's from the inner galaxy and 3C273 above 
atmospheric background for a DUMAND type detector are discussed in the 
cc .text of the Weinherg-Salam model with sin^Qy- 0.2 and Including the 
atmospheric background from the decay of charmed mesons. Constraints on 
cosmological high energy neutrino production models are also discussed. 

It appears that important high energy neutrino astronomy may be possible 
with DUMAND, but very long observing times are required. 

1. Introduction 


There have been a number of recent papers estimating high-energy 
neutrino fluxes and spectra from various astropnysical processes 1 * . 

In this work, I reexamine the problem by 1) pr-'3enting the results of a 
detailed calculation of galactic v-prodvction in cosmic-ray interactions 
with interstellar gas and extragalactlc v-production by interactions of 
ultrahigh energy cosmic rays with the 3K universal background radiation, 
2) using these results t-gether witn present Y -ray observations to pre- 
dict neutrino fluxes and event rates, and 3) discussing the problem of 
observing extraterrestrial v fluxes above the atmospheric background. 

2. Production Rates 


The first basic production process for high energy cosmic neutrinos 
is the decay of charged pions produced in cosmic ray interactions with 
interstellar gns, primarily pp interactions. This process will hence- 
forth be referred to as "pp", although the effects of c/p pHe and orHe 
interactions are included in the calculation. The second process in- 
volves the photoproduction of rr-mesons by interaction of ultrahigh 
energy cosmic rays with the 3K universal microwave background radiation 
(henceforth referred to as yp) and subsequent meson decay. Both of these 
types of interactions involve the accompanying production of tt° mesons 
and their decay into cosmic Y -rays. Thus, the production rates of cosmic 
Y-rays and neutrinos are generically linked. 

2 . 1 pp Neutrino Production 

A detailed discussion of the kinematics of the production and decay 
of secondary particles produced in "pp" interactions may be found in 
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Reference 7. Details of the yp process have also been previously. 

given. ** ) The *-ray production soectrum in v p interactions has been cal- 
culated by several workers. a> 3 > :i * 


In the present paper, the "pp" neutrino production spectrum was 
calculated for pp interactions up to 30 GeV by methods previously 
employed. 7 • 11 * At higher energies, scaling was assumed to hold. At 
these energies, the parameters of Ganguli and Sreekantan 12 * were adopted 
for the rapidity distribution of charged pions. These authors have 
calculated t- ray spectra which are in good agreement with those of the 
present author 13 * for * - ray energies ^ 10 GeV. At higher energies, the 
assumption of scaling gives larger fluxes and a flatter spectrum than 
those calculated using an "isobar + fireball" model in "pre-scaling 
days". 7 * It should be noted that the isobar + fireball (I-F) model is 
equivalent to a "leading-pion" model at high energies, since the isobar 
carries off ~ 50% of the energy and decays into a "leading pion". Both 
scaling and the I-F models produce a secondary spectrum which has the 
same spectral index as the primary spectrum in the high-energy limit. 

The results of tne "pp" production spectrum calculation are shown 
in Fig. 1. The two neutrino production spectra are given for an Inter- 


Fig. 1 - Differential produc- 
tion spectra of neutrinos and 
T-rays from the decay of pions 
produced by interactions of 
cosmic-rays in our galactic 
neighborhood with interstellar 
gas having a mean hydrogen 
density of 1 atom per cm 3 . 

The Y-ray curve and the upper 
neutrino curve are calculated 
for cosmic rays having a spec- 
tral index of 2.67 between 10 
and 3x10 s GeV; the lower neu- 
trino curve is for a cosmic 
ray spectrum with index 2.75. 
The spread in the curve is 
indicative of the uncertainty 
in such calculations. 



LOG E (Gev) 
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stellar hydrogen density of lcm“ 4 so that it is really a production rate 
per hydrogen atom. The upper neutrino curve and the 7- ray production 
spectrum shown in the figure are calculated usin$ a primary cosmic ray 

spectrum Ip(E) ■ 2.35 E~' * 6 7 cm _i s“ l 8r _1 GeV~ 1 for Ep > 10 GeV and the 
lower curve is obtained for a primary spectrum Ip(Ep) ■ 2.0 E“ 2 * ?5 cm” 7 
s -1 sr" l CeV‘ 1 (Ref. 14). (Most recently, Goodman, et al. u 'have reported 
a primary proton spectrum of 1.5 E" ? ’ 7 l *° • 0 *cm -2 s" 1 sr“ 1 CeV“ 1 for l(f>E>l(f 
GeV, consistent with the lower value given In Ref. 14.) The y-ray spec- 
trum is for all y-rays from TT°-decay; the v-spectrum Is for the vp compo- 
nent. Since each pion decay results in one Vy, one v’ u and one V e (V 2 ), all 
of roughly the same energy, the total number of v’s produced (although 
of four different types) is a factor of 3 higher than that shown for \>y'8. 
in this regard, it should be kept In mind that the cross section ratios 
given in Table 1 apply. 


Table 1. yN Cross Section Ratios 

Ratio Experimental Value (Refs. 54,55) Theoretical Value 



1.26 +0.23 

°v N /a V N 

1.32 +0. 32 

e 


a V /0v e N 

0.40 +0.12 

a v e N/ a y e N 

0.38 +0.02 


*valence quarks dominant 


Figure 2 shows the integral 
galactic "pp" neutrino produc- 
tion spectrum. 


Fig. 2 - The integral v-produc- 
tion spectrum obtained from the 
upper neutrino curve in Fig. 1. 
The high energy cosmic ray 
spectrum used is from a recent 
analysis of Hillas (private 
communication) . 


2 . 2 p* Neutrino Product ion 

The production rate for yp 
interactions has been calculated 
using the method of Ref. 9. As 
an example of the range of 
uncertainty, two different 
primary spectra were used, i.e. 

I p (F. ) = 6.4xlO" 31 (E /10 11 )" 2 * 6 
cm" il s _x sr“ 1 GeV -1 (Ref. 16 ) and 


1 

1 

~ 1/3* 
~ 1/3* 



0 
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Ip(Ep) - 2. 4xlO“ Jl (Ep/lO* 1 ) ~ 3 -li *cm" ii s" 1 sr” 1 CeV“ 1 (Ref . 17 ). The latest 
analysis of the ulcrahigh energy air-shower data (Hillaa, private comcn'- 
nicatior.) gives I p (E p ) =- 4xl0“ 31 (E p /10 x 1 ) ■ 3 cm _a s" 1 sr" x GeV _i in the energy 
range 10 a ^Ep6,ll> 0 GeV with indications of a flattening to I p (E ) »* 

2xl0“ 3 ° (Ep/10 11 )" a * 3 in the energy range 1 GeV. The P results of 

the Yp calculation are shown in Figure 3 for the two spectra chosen. The 
right-hand scale of Figure 3 also shows the diffuse background flux from 
this process obtained by multiplying by the factoi c/(4nH 0 ) where H 0 , the 
Hubble constant, is taken to be 50km s~ 1 Mpc“ 1 . This flux, which is also 
shown in the subsequent figures, only holds if the ultrahigh energy 
primary spectrum is universal, an assumption which is contradicted by the 
lack of an observed high erergy cutoff in the spectrum 8 * ' *" 20 1 . Thus the 
fluxes given may be overestimated and may actually be upper limits. 


3. Diffuse Cosmic Neutrino Fluxes 


Figures 4 and 5 show the differential and Integral v^+v^ fluxes 
from various sources (the corresponding flux of electron-neutrinos end 
antineutrinos is lower by a factor of two.) In Figure 4, the cross- 
hatched region marked yp is obtained from the curves shown in Figure 3. 
The hatched region marked pp(G. C. |b |sl0°) is for galactic v-production 
coming from the galactic central region defined by galactic longitude 
330°s^£40° and latitude |b|il0°. The dashed line is the flux computed 
using the isobar-fireball (I-F) model 7 *, which is probably too low. The 
lines bounding the hatched region are obtained using the two "pp" 
production spectra shown in Figure 1. All galactic "pp" curves are 
normalized by assuming that ^ of the galactic y -radiation above 100 MeV 
observed by SAS-2 2 * is from n^-decay and by relating the y-ray and 
v-production using the results given in Figure 1. Spectral measurements 
indicate that ~ 50% of the galactic y-rays observed above 100 MeV are 
most probably from cosmic-ray electron bremsstrahlung with no associated 
v-production. The curve marked PCR (UL) is an upper limit on the v-flux 
from prinordial (or "cosmological") cosmic-ray interactions obtained 
under the assumption that the extragalacti c y-ray background (8xl0' a 
cm” a s“ 1 sr“ 1 above 100 MeV with an E -3 power law spectrum) is from inter- 
actions of cosmic-ravs at high redshifts . 23 * A detailed comparison of 
the calculated and observed spectra 2 ! * indicates that the y-ray back- 
ground is most probably not. from PCRs so that the curve shown is an 
upper limit. Note that this upper limit is many orders of magnitude 
below the absolute cosmological neutrino upper limit given in Reference 
3 (see Section 7) . 

The integral fluxes in Figure 5 show galactic "pp" muon-neutrino 
and antineutrino fluxes for the central region 330°£-(£40 o (GC) for 
|b|^2° and | b | =*10° , the anticenter region (AC) for jb|s2° and the region 
of the galactic poles (GP) . All of these spectra were normalized using 
the SAS-2 y-ray data"'* except for the CP curve which was obtained from 
Fig. 2 using a mean galactic path length of hydrogen gas <nL> = 
3xlO ao cm‘\ The atmospheric Vp+v p flux shown in Fig. 5 is from Allkofer 
et al. 2G ) who used the atmospheric muon data to derive their fluxes. 

The curve marked "other galaxies" is obtained by assuming a density of 
normal galaxies of 3xl0 _7e cm“ 3 (Ref. 27), assuming that they have the 
same y-ray luminosity as our galaxy which is assumed to radiate of its 
Y-ray flux from TT°-decay. Note that since the predicted Y-ray spectrum 
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from normal galaxies Is 
considerably flatter 

than the observed extra- 
galactic background, so 
that the limit on the 
Y-ray flux from other 
galaxies due to n°-decay 
is probably ~ 2x10“® 
cm -1 's“ 1 sr” 1 above 100 MeV, 
the v-background from 
n -decay cannot be sig- 
nificantly higher than 
the curve marked "other 
galaxies". 

Figure 6 shows the 
latitude distribution of 
Vp+V^ calculated for 
energies above 1 TeV by 
extrapolation from the 
SAS-2 y-ray data. 



Fig. 3 - Calculated differential neutrino 
production spectra and background fluxes from 
photopion production by ultrahigh energy cosmic 
rays interacting with the 3K universal micro- 
wave background radiation. 


4 . Prompt Atmospheric Neutrinos 


A new source of "prompt" neutrinos, possibly from the decay of 
charmed mesons, has been fotmd in recent beam dump experiments at 
CERN 26 “ ,0) . The ratio of the production rates of these neutrinos to 
TT-production is R^IO 4 . A summary of the present data on charmed meson 
production is given in Reference 31. The importance of prompt neutrino 
production in the atmosphere at high energies has been discussed by 
several workers. 3 * * 2, 3S ' 


Figures 7 and 8 show the estimated fluxes of prompt (PR) muon and 
electron neutrinos for R •- 10“ 4 and R - 10 -3 (an R value between these 
two numbers appears most likely at this writing) together with the 
horizontal and vertical fluxes of atmospheric neutrino3 2?J and our 
estimates of the neutrino flux from the central part of the Galaxy. 

It can be seen that prompt neutrinos from the atmosphere have a spectrum 
and intensity which may mimic those of galactic neutrinos. Furthermore, 
even with directional information, a "large" background of prompt 
neutrinos may make it difficult to pick out a galactic "signal". 

The neutrino flux from the inner Galaxy is expected to be confined 
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mainly to a narrow band in the sky (see Fig. 6) which, in principle, can 
be picked out above the atmospheric background’*' with a detector such 
as the proposed DUMAND system with an angular resolution of better than 
1/2° above 10 3 GeV energy. 1 (see section 6.1). 

5. Event Rates 


The vN cross section rises linearly with E v and is ~ 0. 7xlO" 36 E v cn^ 
at accelerator energies. The vM cross sections are ~ 1/3-0. 4 of this 
value (see Table 1), although there is now evidence of an increase in 
this ratio to ~ 0.6 at E y ~ 100 CeV ,7, t possibly due to a new flavor fa- 
quark with mb ~ 5 GeV.’** Owing to scaling violations and quantum 
chromodynamics, this ratio may approach unity at E^ 10 6 GeV. s ** Above 
l critical en« rgy 

E V C ~ Mw 2 /2Mp. (1) 


where My is the mass of the intermediate vector boson and Mp is the 
proton mass, the energy dependence of the cross section levels off to a 
logarithmic one. 40) In the unified gauge theory of Weinb :g and Salam, 
My ■ 37.3 esc 8y, where 9y is the mixing angle parameter expressing the 
relationship between the neutral current and electromagnetic coupling 
constants. 4 1 The best experimentally determined value at present is 
sin a 9y =>» 0.2 (Ref. 42) giving My “• 84 GeV. The resulting vN cross sec- 
tion as a function of E y is shown in Fig. 9. 


Using the vN cross 
section from Fig. 6 for My * 
84 GeV (solid line) and a 
strictly linear cross section 
(My - 00 dashed line) together 
with the calculated cosmic 
v-f luxes, expected event 
rates for v 's falling on a 
detector with an effective 
mass of 10° tons of seawater 
as proposed for the DUMAND 
experiment were calculated 
and are given in Fig. 10. 

The dotted line narks the 
level corresponding to 100 
events/yr. The upper "pp" 
curve is for the flux from 
the galactic central region 
using the production spec- 
trum of Fig. 2, which is 
derived from the upper v- 
curve of Fig. 1. The curve 
marked <PP>(-al * s ^ or t * ie 
galactic "pp ^neutrinos 
averaged over 4n sr. Thus, 
on the average, one would 
expect an event rate of ~ 40 
yr -1 from galactic v^'s of 
energy > 10 TeV. If ultra- 



Fig. 4 - Differential background v-fluxes. 
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Fig. 5 - Integral Vp+Vp fluxes 
from cosmic sources and atmos- 
pheric tt and K decay. 



Fig. 6 - The galactic latitude 
distribution of neutrinos with 
energy above 1 TeV in the direc- 
tion of the inner galaxy. 


high energy cosmic rays are universal (not likely) and the vN cross 
section rises linearly up to ~ 10 9 GeV (My^5xlO*GeV, again not likely) 
this rate could again be reached for E v > 10 9 GeV from Yp interactions 
except for the fact that if na v ^ becomes greater than the effective 
geometric area of the detector, A, the detector becomes area limited and 
the event rate is IyA s -l sr _1 (see Fig. 10). (However, in this case, the 
ocean above a DUMAND type detector would be opaque to v's.) It is more 
likely, however, that at energies high enough to see sceady diffuse 
fluxes of cosmic neutrinos with a 1(7° ton detector above the atmospheric 
background, the event rates will be very low. The resulting statistics 
will make it hard to see galactic neutrinos above an isotropic back- 
ground of prompt neutrinos (see section 6.1). 


6 . Observability of Extraterrestrial Sourc es 

The large fluxes of high energy neutrinos produced In the atmosphere 
as the result of the decay of secondary mesons arising ir. cosmic-ray 
interactions with atmospheric nuclei provide a significant background 
which, although useful from the point of view of neutrino physics, pro- 
vides constraints on conducting neutrino astronomy observations. Below 
~ 10 TeV, neutrinos from tt and K decay create a noisy background problem 
which, because of the softness of their spectrum, gets worse at lower 
energies. At energies above ~ 10 TeV, the prompt neutrinos from D-decay 
take one more significance because of their harder spectrum which mimics 
the galactic spectrum (see section 4). One must therefore consider the 
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observability of various cosmic neutrino sources above fluctuations in 

this noisy atmospheric background. 


6.1 Observability 
Galaxy 


of the Inner 


The observability criterion 
adopted for quantitative estimates 
is the number of standard devia- 
tions a cosmic source signal 
would provide above atmospheric 
fluctuations. Based on the 
results given in the previous 
section, we find that the number 
of galactic neutrinos (v^+v^) 
observed in «.ne year at E i 10 TeV 
in a 10° ton DUMAND detector is 
calculated to be 45 if one used 
a primary cosmic ray spectrum 
dJ/dE “ E” 3 ' 67 . However, adopting 
the steeper primary spectrum 
obtained in more recent work 1 ** 16 
a more conservative event rate is 
obtained. The primary spectrum 
given in Ref. 14, which is consis- 
tent with the data in Ref.’s 15 
and 26, will be adopted here, 
yielding a galactic event rate of 
~ 18/yr above 10 TeV. At 1 TeV, 
the event rate is ~ 16 times 
higher. 



F.. (G«v) 


About 45% of the expected 
galactic events come from a 
galactic central region +8° in 
latitude and +50° in longitude. 
»*he prompt atmospheric neutrino 
fluxes for R * 10“ 4 and R “ 10“ 3 


Fig. 7 - Fluxes of atmospheric muon- 
neutrinos from tt and K decay and from 
prompt decay as compared with the 
predicted diffuse flux from the inner 
galaxy. 


are given in Figure 7. The atmospheric n and K decay fluxes are given 


in Ref. 26. The mean atmospheric neutrino flux from n and K decay is 
~ 0.5 of the horizontal flux at 1 TeV and ~ 0.3 at 10 TeV.‘ ,3j 


The event rates for a 10° ton detector as given in Ref. 34 are shown 
in Tablt 2. It can be seen that the inner region of the galaxy may be 
detectable above atmospheric background fluctuations with observing runs 
of 4 years or longer. 


6.2 Observability of the Nearby Quasar 3C273 

In considering neutrino astronomy of extragalactic point sources, 
the quasar 3C273 provides an excellent subject for consideration because 
of its closeress and prominance among quasars and because it is the only 
knovn extragalactic source of 100 MeV v-rays. Thus, if the Y-radiation 
from 3C273 is the result of meson production in cosmic-ray interactions 
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Table 2. Observability of Galactic Central Region 



Energy 

N_ 

Prompt 

(R-10 -3 ) 

380 

Number 

of Atmospheric 

Number of 

c's of 

V. +v 

(TeV) 
> 1 

Neutrino Events/yr. 
»vr. 

17.600 130 

v's from 
G.C./yr . 

130 

G.C. Events 

Above Atm. 
Bkgd. (4 yr) 

2.0 

V V 

>10 

23 

120 

12 

8 

1.4 

v e +v e 

> 1 

380 

790 

34 

65 

3.8 


>10 

23 

4 

5.2 

4 

1.5 

V 

> 1 

(R-10 -4 ) 

38 

17,600 

130 

130 

2.0 

V U 

>10 

2.3 

120 

11 

8 

1.4 


> 1 

38 

790 

29 

65 

4.5 


>10 

2.3 

4 

2.5 

4 

3.2 


In the source, either through tt°- 
decay or the breysstrahlung of 
electrons from it decay, the 100 
MeV neutrino flux from 3C273 may 
be comparable to the ^6x10 7 
cm 2 9 ' y-ray flux observed by 
the COS-B satftllite 44, for this 
source. 

Taking the resolution element 
for neutrinos with E>1 TeV from 
DUMAND to beSl R - 3xl0' 4 sr ,f) , the 
atmospheric neutrino event rate In 
is expected to be 'vll/yr above 
1 TeV. Thus, a 3o signal, certainly 
the absolute minimum for observa- 
bility of a point source, would be 
'''lO/yr. For a source with a 
differential energy spectrum 
“ E 2 * 7 , this would correspond to 
a il_TeV neutrl to f lux of — 6xl0“ ia 
cm 2 s l . If the integral v-spectrum 
from 3C273 follows an F ^‘ 7 gower law 
from a value of 6x10 ’cn 2 s 1 
at 100 MeV, the resultant flux above 
1 TeV would be '‘-10 12 co 2 s 1 , a 
factor of 60 bel w the level of 
detectability. Indeed, the Integral 
spectrum would have to be flatter than 
E 1,25 for detectability, which seems 
unlikely since the photon spectrum 
appears to steepen between the x-ray 
and y ray range 4 ®. The case for 
detectability of other steady extra- 
galactic sources whose 100 MeV 



Fig. 8 - Fluxes of atmospheric 
electron-neutrinos from tt and K 
decay and from prompt decay as com- 
pared with the diffuse flux from 
the inner galaxy. 
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y-ray flux remains below present 
satellite detectability seems 
even Tore hopeless with one 
possible loophole. It is possible 
that cosmic rays could be produced 
in the core regions of quasars, 
there producing v's while the 
Y-rays are absorbed in optically 
thick regions. In that case, the 
producing regions must be much 
smaller than the <0.1 pc size 
of the x-ray producing region 
since the column density in the 
x-ray producing regions is 
< 4.5xl0 ai cnF at the 90% confi- 
dence level * less than 4 orders 
of magnitude below the Y-ray meart- 
free-path. It also appears 
that the Y-ray emission from 
quasars is more likely g elec tro- 
magnetlc in origin 47 ’ 9) and 
therefore unconnec ed to v-pro- 
duction. It would therefore 
be of extraordinary interest if 
neutrinos from 3C273 were 
detected, as it would reveal 
entirely new information aoout 
the fundamental nature of quasars* 
But observability is uncertain 
at best in terms of present theo- 
retical ideas. 



Fig. 9 - Cross section for vN inter- 
actions assuming an intermediate vector 
boson mass My ■ 84 GeV. 


7. Cosmological Neutrino Background 

Following a suggestion 49 * that a high-redahift cosmological origin 
of the high energy cosmic radiation could explain the steepening of its 
spectrum above ~ 10“ GeV, a cascade origin of the cosmic Y-ray background 
radiation was also suggested 50 *, followed by suggestions that such a burst 
of cosmic ray production in the distant past would result in copious 
neutrino fluxes 9 It should, however, be noted that photopion produc- 
tion of neutrinos in an early intergalactic (or pregalactic) medium is 
limited by observational constraints on the associated Y-ray background 
radiation 7 *. These constraints therefore call into question burst pro- 
duction or "bright phase" models. The "pp" limit from the Y-ray obser- 
vations is shown in Fig. 5 marked PCR(ITL) (see section 3). The py limit 
is not related to the x-ray background av previously suggested 51 * because 
the x-ray background cannot be related to photopion production since the 
shape of the Y-ray background has the form ®F. -3 above 10 MeV ?4 * whereas 
the photopion model would give a spectrum of the form ' 0> “F." a . Further- 
more the model on which the high-redshif t calculations are based Implies 
a Y-ray background flux above that observed unless the mean intergalactic 
gas density at present is < 10" ,0 cm -3 (based on the discussion in Ref. 23 
with updated y-ray data) , and unless galaxy formation occurred at an 
extremely efficient rate at redshifts zsl5. These assumptions create 
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problems of compatibility with 
recent x-ray observations of inter- 

galactic gas in galaxy clusters 
and with radio and optical galaiqr 
observations 5 2) . 

The remaining way of produc- 
ing a high flux of neutrinos at 
high redshlfcs would be In the 
cores of compact astrophyslcal 
objects (quasars, Seyfert 
galaxies, etc. *')• Compatibility 
with other observations would be 
maintained provided that a) the 
co-produced Y-radiatlon is 
absorbed In the sources, and b) 
the cosmic rays also do not escape 
the sources and produce lnter- 
galactlc v-ray emission. The 
difficulties with such models 
have already been discussed In 
section 6.2. 
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